A consequence of the gravitational self-force for 
circular orbits of the Schwarzschild geometry 
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A small mass ^ in orbit about a much more massive black hole m moves along a world line that 
deviates from a geodesic of the black hole geometry by 0{^/m). This deviation is said to be caused 
by the gravitational self-force of the metric perturbation hab from /i. For circular orbits about a 
non-rotating black hole we numerically calculate the 0{fi/m) effects upon the orbital frequency and 
upon the rate of passage of proper time on the worldline. These two effects are independent of the 
choice of gauge for hat and are observable in principle. For distant orbits, our numerical results 
agree with a post- Newtonian analysis including terms of order {v/cf . 

PACS numbers: 04.25Nx, 04.30.-w, 04.80.Nn, 97.60.Jd 



I. INTRODUCTION AND SUMMARY 



Gravitational waveforms, perturbation analysis 
and the self-force 



The push towards the detection of gravitational waves 
depends upon accurate theoretical models of the sources. 
Numerical relativity has made great progress in the past 
few years and now tracks black hole binary systems for 
many orbits and down to the final coalescence and ring 
down. These recent computational successes arc joining 
the post-Newtonian analyses in providing building blocks 
for the construction of gravitational wave templates. 

A gap in the theoretical progress remains for binary 
systems with an extreme mass ratio. A small black hole 
in a close orbit about a much larger one might be too 
relativistic for post-Newtonian analyses. And the great 
difference in length scales for the two black holes causes 
difhculty for a full numerical relativity solution. This gap 
is the natural domain of perturbation theory. 

The state of the art in pcrturbative gravitational wave- 
forms [1] has a small mass /i moving along a geodesic of 
the metric gab of the much larger black hole. The per- 
turbed Einstein equations then determine the metric per- 
turbation hab which contains the waveform. But, these 
waveforms would be much improved if the influence of 
hab back on the motion of were part of the analysis. 
Then the dissipative effects of radiation reaction, joining 
with conservative effects, would be reflected in the wave- 
form and would allow tracking for more orbits without 
losing phase information as the inspiral proceeds. 

The metric perturbation hab causes the worldline of /i 
to be a geodesic of gab + hab- Sometimes the same motion 
is described as non-geodesic in gab with the acceleration 
Cgp being caused by the self-force from hab acting on 
/i. In this paper, we prefer the "geodesic in gab + hab" 



description of the motion, and use the phrase "self-force" 
only in a generic way to describe any and all of the 0{fi) 
effects of hab on the world-line of fj,, even those which are 
unrelated to the acceleration Ogp. 

When the mass ^ is modeled as a point particle, an 
attempt to calculate any consequence of the self-force 
immediately bumps against a technical difficulty. The 
metric perturbation hab is singular at the location of ^, 
precisely where it must be evaluated. Some method of 
regularization must be invoked to remove the singular 
feature. The formal aspects of this regularization are now 
very well understood [2-13], and we give some details for 
our particular application in Sec. IV. 



B. Our limitations 

In this paper we limit our analysis to the self-force 
effects upon the circular orbits of the Schwarzschild ge- 
ometry. The numerical analysis is straightforward, and 
the theoretical interpretation of the results is clear and 
unambiguous. These simple orbits provide a convenient 
test bed for ideas regarding the self-force. 

Despite the emphasis of our opening paragraphs, this 
paper contains no new information regarding the dissi- 
pative effects of radiation reaction on the gravitational 
waveform for slow, circular inspiral into a non-rotating 
black hole. In Appendix D we sketch a proof, based 
upon the form of the perturbed Einstein equation, that 
the self-force removes energy directly from the orbit of 
11, via Eq. (19), at a rate that matches the energy be- 
ing radiated away in the wave zone and down through 
the event horizon. This result is satisfying but not unex- 
pected. At this order of approximation, the dissipative 
effects of the self-force imply that the orbital frequency 

should change at a rate 
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where dEjdt is evaluated from Eq. (19), and y,E is the 
0(/i) contribution to the mass-energy of the system as 
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measured at a great distance, which is equivalent to 
fiT''Ua + 0{fj.'^) \ with being the timeUke Kihing vec- 
tor and being the four velocity of fi. Thus, for circular 
orbits of the Schwarzschild geometry the self-force for- 
malism based upon first-order perturbation theory adds 
formality but yields nothing new about the evolution of 
the gravitational waveforms from the dissipative part of 
the self-force. In fact, conservative self-force effects on 
the values of E and for circular orbits are described 
in Sec. II B but these are not sufficient for improving the 
calculation of cKl/dt from dissipation without at least 
also determining the 0{^^) part of dE/dt. This result is 
common knowledge within the self-force community [16] 
but appears not to have spread much beyond that. 

Non-circular orbits in the Schwarzschild metric and all 
orbits in the Kerr metric contain substantial complica- 
tions that we also do not address. 



C. Our results 
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FIG. 1: The quantity iit*, which is the gauge independent 
0(/i) part of w* effected by the self- force, is given as a function 
of Rn for circular orbits in the Schwarzschild geometry. Also 
shown are iit* as calculated with IPN and 2PN analyses of 
Appendix E based upon Refs. [17] and [18]. 
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In perturbative matters, we prefer to focus upon those 
observable quantities which can be defined in a manner 
independent of the gauge choice for hab- The angular 
velocity Q and the time component of the four velocity 
M* are two such quantities for /i in a circular orbit. Our 
main computational result is determining the change in 
the relationship between U and u* caused by the self- 
force. The coordinate radius of the orbit is not included 
in this analysis because it depends upon the gauge choice 
and has no inherent physical meaning. 

First, we define Rfi via il^ = m/Rf^ as a natural ra- 
dial measure of the orbit which inherits the property of 
gauge independence from fl. In Sec. Ill we show that it* 
naturally separates into two parts it* = qu* + iit*, where 
each part is individually gauge independent. Further the 
functional relationships between gw* and i?j2 are iden- 
tical to their relationships in the geodesic hmit, fi ^ 0, piG. 2: The same as Fig. 1, but including the 3PN analysis, 
and show no effect from the self- force. 

The remainder im* is, however, a true consequence of 
the self- force, and we plot lU* as a function of i?si in 
Figs. 1 and 2. These two figures along with Table I are 
the primary numerical results of this paper. 

In Appendix E we derive a post-Newtonian (PN) ex- 
pansion for lit* based upon the work of others [17, 18]. 
Our expansion is in powers of m/Rn, which is j(? in 
the Newtonian limit, and from Eq. (E26) 



1 1 1 1 1 1 1 1 1 
• 


1 


1 










1 PN 






-- 2PN 




• 


— 3 PN 




— \ • 






























1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



10 15 

Ra/m 



20 



m 



m 



Rn 



Our Figs. 1 and 2 display the different levels of the 
post-Newtonian expansion. These agree with our numer- 
ical results as well as they should. 

With a good bit of hubris, we have gone one step 
further. After the three known post-Newtonian terms 
were removed from our determined values of im*, we nu- 
merically fit the residual to an expansion of higher or- 
der terms, and found the coefficient of (m/i?si)^ to be 
-27.61 ± .03 2. 

We now describe an astronomical thought-observation. 



which includes terms of order j 



(2) 



^ This equivalency follows from Zerilli's [14, 15] analysis of 
monopole perturbations of the Schwarzschild metric. 



^ We had already used a similar technique while only knowing the 
analytic expansion up to the (m/i?n)^ term. Finding the coef- 
ficient of {m/Rci)^ to be —5.001, with the last digit uncertain, 
provided the impetus to push the analytical work to one more 
order. 
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Imagine that an observed system has a small mass in a 
circular orbit, with orbital frequency f2, about a large 
non-rotating black hole. Assume that the mass m of 
the black hole could be determined by some independent 
means, and that a redshift observation, perhaps as de- 
scribed in Sec. II C 2, could determine u*. Then 

ou' ^[l~i{nmf/^]-^/^ (3) 

from Eq. (42), 

Rn - (^) (4) 

from the definition of i?si, and 

— — qu*. (5) 

With values of m, Rn and lU* in hand, an astronomer 
could then determine the mass of the small object from 
the graphical data. 

Our self-force result may only be of academic interest. 
But it is a result. Our numerical calculations are con- 
sistent with a post-Newtonian expansion. We estimate 
the unknown coefficient of v^/c^ in the expansion. And 
we have proposed one unlikely observational method to 
infer the mass of a small object which is orbiting a much 
larger black hole. 

D. Outline 

In Sec. II we introduce our notation, and we provide 
the components of the geodesic equation for the per- 
turbed Schwarzschild metric gab + hab in Sec. II A. 

In Sec. II B we focus on quasi-circular orbits — those 
that would be circular except for the dissipative effects 
of the self-force. And we find that two components of the 
four- velocity, u* and u^, are effected by the self- force in 
a manner that is 0{fi) and independent of the particular 
gauge that is chosen for hab- Their ratio is the orbital 
frequency il = u'^/u^, which is also gauge independent. 

The component m* is the ratio of the rate of passage of 
Schwarzschild time to proper time at the particle. Two 
different physical interpretations of u* for quasi-circular 
orbits arc given in Sec. II C; one of these is in terms of a 
redshift measurement, the other is related to the helical 
Killing vector field of the perturbed geometry. 

The formal effect of the self-force on the relationship 
between f2 and is described in Sec. III. Particular 
emphasis is placed upon how the gauge independence of 
this effect reveals itself. 

Section IV contains details of our numerical determi- 
nation of the relationship between m* and fl. Table I 
contains an extensive summary of our numerical results 
which might be used by others for a numerical compari- 
son. This section also has a description of the regulariza- 
tion required by our application which points out some 
difficulties inherent in regularizing any quantity which is 



not gauge independent, such as the acceleration of fi as 
measured in the background metric. 

Our conclusion in Sec. V describes both the current 
state of self-force calculations as well as its possible fu- 
ture. 

We have relegated a number of side discussions to Ap- 
pendices. 

Some convenient identities are given in Appendix A. 

Gauge transformations that lead to the understanding 
that u*, u*^ and f2 are gauge independent are given in 
Appendix B. 

An elementary, completely Newtonian example of the 
confusion that gauge freedom can cause is in Appendix C. 

In Appendix D we show that the rate at which /x loses 
energy via the self-force precisely matches the rate that 
gravitational waves carry energy out at infinity and down 
the black hole. 

The post-Newtonian expansion for in Eq. (2) is 
derived in Appendix E. 



II. GEODESICS OF THE PERTURBED 
SCHWARZSCHILD GEOMETRY 

Self-force analysis implies that a point mass fi moves 
along a geodesic of the perturbed Schwarzschild metric 
dab + hab where, in this section, hab is the suitably reg- 
ularized metric perturbation and is at least C^. We as- 
sume that the geodesic is in the equatorial plane. We 
assume that hab = 0{ii), and we work only through first 
order in ji or hab- Indices are lowered and raised with 
gab + hab and g""^ — h""^. We use Schwarzschild coordi- 
nates, and hab is a function of t, r, 9 and 0. To be clear: 
we specifically choose to normalize the four velocity of 
the particle with respect to the actual regularized physi- 
cal metric gab + hab, and not with respect to the abstract 
background Schwarzschild metric gab- 

The quantities E, R, and J are functions of the proper 
time s and are defined in terms of the particle's four- 
velocity by 

E = —Ut, R = and J = u^, (6) 

which leads to 

/ R + u^hb'' \ 

and 

\ 1 — 2m/r / 

where the overdot represents d/ds. Thus, u°' depends 
upon s implicitly through its dependence upon E, R, J 
and also upon r, evaluated at the particle. These expres- 
sions for Ua and u" are recursive, but the recursion is 
intended to be carried only through first order in hab- 
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It is convenient to define the non-radial part of the 
four-velocity of the particle, with no hab terms included, 
as 



E 



J 

,0,0,- 



(9) 



1 — 2m /r ' ' ' 
Also, the instantaneous angular velocity of the particle is 

n = u'f'/u\ (10) 



result is 

R{1 + h/) + -R^drh/ + u'^vl'dahb' + Ehtr + i/i'^'' 



--(l-u'^i2X6) + ^(r-3m) 
+ (l-2m/r)ia.(ruXb). (16) 



B. Quasi-circular orbits 



A. Geodesic equation 

One convenient exact form of the geodesic equation for 
the metric gab + hab is 



dua 1 b c ^ / , . N 



(11) 



While analyzing this geodesic equation we drop all 
terms of 0{(j?). The components of Eq. (11) provide 



ds 



2 dt ' 



d / rR 



dJ_ _ I a bdhab 
ds ~ 2 dcj)' 
d 



(12) 
(13) 



ds\r — 2m 



\ 1 d 

+ U''har) = T;U°'u'' — {gab + hab).iU) 

/ 2 or 



These and other equations describing the particle's 
worldline are to be evaluated at the location of the par- 
ticle. 

The normalization of m° is a first integral of the 
geodesic equation. 



1 = -u°-u'{gab + hab) 



'— + U''v!'hab - R'-hr 



(15) 



1 — 2m I r 

which is derived by beginning with Eq. (8). 

1. The radial component 

Further analysis of the radial component of the 
geodesic equation is lengthy and requires care. We begin 
with Eq. (14), and use Eqs. (7) and (8) extensively along 
with the elementary identities in Appendix A which allow 
for the easy commutation of ?2° with dr- Simplification 
proceeds by identifying all terms which are 0(/i^) if the 
motion is quasi-circular as is described in Sec. II B; these 
terms are systematically moved to the left hand side. To- 
ward the end of the analysis the normalization condition 
(15) is used to replace E'^ in favor of on the right 
hand side. Final simplification arises from a convenient 
placement of the indices of hab on the left hand side. The 



For any orbit Eqs. (12) and (13) imply that E ^ J ^ 
0(/i). In this section our interest focuses on the orbits 
which would be circular but for these effects of radiation 
reaction. For these quasi- circular orbits we require that 
R = 0{^) and also that R = 0{iJ?). The second of these 
conditions is necessary to insure that the orbit is not a 
precessing elliptical orbit with a very small eccentricity 
e = 0{fi). To summarize, for a quasi-circular orbit E, R 
and J change slowly, E ^ R ^ J ^ 0{fi), and their rates 
of change are slower still, E ^ R ^ J ^ 0{iJ?). 

In this section wc continue to drop 0{ii^) terms includ- 
ing those which arc consequences of the quasi-circular 
condition from the geodesic equations. Thus the normal- 
ization of becomes 



1 = 



E' 



J' 



1 — 2m I r 



-f u°-u^hab. 



(17) 



Neither T^d/dx"- = d/dt nor <i>''d/dx'' = 9/90 is a 
Killing vector of gab + hab, but a helical Killing vector fc" 
exists for a quasi-circular orbit in that Ck{gab + hab) = 
0(/i^). The Killing vector exists for any gauge choice for 
hab- But if the gauge choice is reasonable and hab re- 
spects the helical symmetry, then is particularly sim- 
ple [13] with 



k^d/dx" = d/dt + nd/dcj), 

Ckgab = 0, and Ckhab = 0{\j?) 



(18) 



In what follows, the only assumption which wc make 
about the gauge is that it respects both the helical sym- 
metry over a dynamical timescale and also the natural 
reflection symmetry through the equatorial plane. Our 
results are independent of the assumption, but the anal- 
ysis is much simplified because of it. For example, ii° is 
tangent to a trajectory of fc" so that u°-dahbc = 0{ij?) at 
the particle for a quasi-circular orbit. 

For quasi-circular orbits Eqs. (12) and (13) simply im- 
ply that 



'^-^ 1 -a-hn . J 

— = -7— tU u dthab and 
at 2u^ 

dJ 1 -a-b^ 1 

- ^ —,UU d.hab, 



and Eq. (18) then yields 

dE 
'dt 



dJ 
'dt' 



(19) 
(20) 

(21) 
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We show in Appendix B that the right hand sides of Eqs. 
(19) and (20) are both independent of the gauge choice 
for hah- Thus, dE/dt and dJ/dt are gauge independent 
for quasi-circular orbits. 

For a quasi-circular orbit every term on the left hand 
side of the radial equation (16) is 0{fj,^). It easily follows 
that 



Further, the ratio of and u* gives the orbital fre- 
quency 



3m 



2r2 



u°-u''drhab (28) 



■11 



-(1 ~ - (^-4^) ^ 

\ \r — 6m/ 2 or ^ ' 



3m ^ \r — 3mJ2dr 

Eqs. (17) and (22) together imply that 



(22) 



E'^r{r - 3m) 
(r — 2m)2 



1 



r d 



u''v'hab---^{u''v'Kb). (23) 



It is important to note that the 0(/i) parts of E and J de- 
pend upon the choice of gauge and have no precise physi- 
cal meaning in terms of the energy or angular momentum 
of the particle or the system. There is some irony that 
dE/ds and dJ/ds are gauge independent, while E and J 
arc not. 

Two interesting quantities for quasi-circular orbits are 
components of the four velocity 

( = f^V = (-E + u'^htb)^ 

' \ds) (1 - 2m/r)2 



3m 



+ - ^U'^u'drhab) (24) 



and 



for a quasi-circular orbit. Its gauge independence is in- 
herited from u* and u'^, and may also be demonstrated 
directly following the same steps as for u*. 

While Q is observable from infinity and is gauge inde- 
pendent, the Schwarzschild radial coordinate of the orbit 
is neither. A gauge transformation simply changes the 
the radius of the orbit by Ar = S,^ . 



C. Physical interpretations of 

The quantity u* is the ratio of the rates of change of 
Schwarzschild coordinate time and of proper time along 
the particle's geodesic in the perturbed and regularized 
Schwarzschild metric. We have two different physical in- 
terpretations which give further meaning to u*. 



1. as a constant of motion 

The perturbed metric gab + hab has a helical Killing 
vector fc". With the reasonable gauge choice that hab 
respects the helical symmetry as in Eq. (18), it follows 
that 



\ ds ' r 
(r - 2m) 
r(r — 3m) 



u^'h^b? 



-k°-ua = E~nj 



(29) 



—— - -U'^u''drhab 

r[r — 2m) 2 



(25) 



is a constant of motion over a dynamical timcscalc. For 
a general orbit 



These follow from Eqs. (23) and (22) along with the equa- 
tions in Appendix A. 

A gauge transformation is an infinitesimal coordinate 



{E - nj)u* = Eu* - Ju'>' 

E{E + u^'htb) J 



transformation : 



■-old 



+$,"', which transforms hab 



hab + ^hab where Ahab = -2V(q4). For hab respecting 
the symmetry of the helical Killing vector, we show in 
Appendix B that A{u'^u''hab) = at /i. Thus the right 
hand side of Eq. (24) changes only in two places. The 



1 — 2m/r T 
= 1 ^ Ru'^har + R" 



— {J - u h^b) 



r — 2m 



.r)30) 



change in r/(r — 3m) resulting from r 

3m 



■C is 



with Eqs. (8) and (15) being used to derive the second 
and third equalities, respectively. For a quasi-circular 
orbit, the right hand side of the last line is 1 + 0{iJ?), 
and thus 



— 3m 
and the change 

Aiu'^u'^drhab) 



(?■ — 3m)- 



6m 

r'^{r — 3m) 



(26) 



(27) 



E-QJ = 1/w* 



(31) 



is from Eq. (B15). These two changes cancel each other at 
0{fi) on the right hand side of Eq. (24), which is therefore 
gauge independent. A similar argument shows that u'^ is 
also gauge independent. 



The gauge independence of then implies that the con- 
stant of the motion E — flJ is also gauge independent 
even while E and J are not. 



2. u' as a redshift measurement 

In principle u* might be measured by a redshift obser- 
vation, which we now describe. 
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Let a light source be located at the particle where the 
perturbed metric is suitably regularized. Let be the 
tangent vector to an afhnely parameterized null geodesic 
of a photon from this light source. The energy f cm of the 
photon, as emitted, is proportional to u'^^Va with it"-,-^ 
the four velocity of the emitter. Thus the ratio of the 
energies as measured when observed and when emitted 
is 



and 



(32) 



with u'^^ the four velocity of the observer. With k°- being 
the helical Killing vector field of the perturbed metric, 
k°-Va is a constant of motion along the geodesic of the 
photon. The four velocity of the emitter is the four veloc- 
ity of the particle, and u'^^ oc k°- so that = u*fc°|em- 
Let the photon be observed along the z-axis at a large dis- 
tance away from the black hole so that u'^^ = 'u-^oo^t = ^t- 
We then have 



fob 



U*-{k°-l'a) 



u^{k°-iyaY 



(33) 



Finally, this may be written as 

gob ^ 



oo 

t 



(34) 



because = at a large distance along the z-axis. 
Thus, the gauge independent determines the red- 

shift of a photon, emitted from the particle, when the 
photon is observed on the z-axis at a large distance. 



III. GAUGE INDEPENDENCE OF OUR 
RESULTS 

In our analyses, we have a gauge independent expres- 
sion for the orbital frequency in Eq. (28), 



1 



r(r - 3m) _a-ba i. 
4to 



(35) 



which allows us to define a related gauge independent 
measure of the separation between m and fi via 



(36) 



iu* = (l-3m/r)-i/2irsX6. 



(39) 



From the discussions above and in Appendix B, it is clear 
that each of these parts of u* are individually gauge in- 
dependent. 

For the moment, assume that the analysis is being 
performed in a special gauge where vf^u^drhab = 0; the 
gauge vector required to put a general hab into this spe- 
cial gauge is given in Eq. (B15). In this case 



f7 



m 



and 



ou*-(l-3m/r)-i/^ 



(40) 



(41) 



Eqs. (40) and (41) are equivalent to their geodesic limits. 
It follows that, for this particular choice of gauge, 



= [l-3(r!j7i)2/3]-i/2 



(42) 



Every bit of this equation is gauge independent. There- 
fore Eq. (42) must hold for any choice of gauge. An al- 
ternative, direct verification of this equation results from 
the use of Eqs. (35)-(38) to write both sides of Eq. (42) 
in terms of r, to, u"' and hab, followed by the expansion 
of each side through linear order in hab- 

From our perspective the relationships between 17, 
and qu* are as if the motion were geodesic in gab and 
independent of the self-force and should be treated as 
such. 

To see the effect of the self- force on w*, we numerically 
determine u'^vJ'hab and, thence, im* as a function of Rn 
as described in the next section. We have no need to 
evaluate vf^u^drhab- 

Note that if the motion of /i were being described as 
"accelerated in the Schwarzschild geometry", then u° 
would be normalized by gab- Also u* would not be gauge 
independent when defined in that manner and the gauge 
independence of m* with our normalization from gab + hab 
would likely be undiscovered. 



IV. DETERMINATION AND 
REGULARIZATION OF hab 



In the geodesic limit, /i — > 0, the separation _Ro is pre- 
cisely the Schwarzschild radial coordinate of the orbit. 

Our other gauge independent quantity of interest is u*, 
in Eq. (24), which we now choose to write as 



The major results or our numerical analysis were dis- 
played in Figs. 1 and 2 of Sec. I C. In this section we fill in 
some details regarding how these results were obtained. 



where 



ou* = (1 - 3m/r)-i/2 _ ''-u'^vl'drhab) 



(37) 



(38) 



The retarded field 



We numerically determine the retarded metric pertur- 
bation h^^^ in the vicinity of a small mass /i in a circular 
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-0.1029412 


1.871470 


X 


10" 


6 


25 


1.066003 


-0.04359984 


-7.98 


X 


10" 


-5 


-0.08180055 


6.157919 


X 


10" 


■7 


30 


1.054092 


-0.03577831 


-3.75 


X 


10" 


-5 


-0.06788457 


2.486475 


X 


10" 


■7 


35 


1.045825 


-0.03034064 


-1.99 


X 


10" 


■5 


-0.05802241 


1.155513 


X 


10" 


-7 


40 


1.039750 


-0.02633967 


-1.15 


X 


10" 


-5 


-0.05066538 


5.950154 


X 


10" 


-8 


50 


1.031421 


-0.02084465 


-4.65 


X 


10" 


-6 


-0.04041928 


1.962457 


X 


10" 


■8 


60 


1.025978 


-0.01724759 


-2.22 


X 


10" 


-6 


-0.03362174 


7.926444 


X 


10" 


■9 


70 


1.022142 


-0.01470964 


-1.19 


X 


10" 


-6 


-0.02878197 


3.681881 


X 


10" 


■9 


80 


1.019294 


-0.01282296 


-6.94 


X 


10" 


-7 


-0.02516046 


1.894535 


X 


10" 


-9 


90 


1.017095 


-0.01136531 


-4.32 


X 


10" 


-7 


-0.02234857 


1.054112 


X 


10" 


■9 


100 


1.015346 


-0.01020528 


-2.82 


X 


10" 


-7 


-0.02010207 


6.238203 


X 


10" 


-10 


120 


1.012739 


-0.008475251 


-1.35 


X 


10" 


-7 


-0.01673728 


2.515767 


X 


10" 


■10 


140 


1.010889 


-0.007246793 


-7.30 


X 


10" 


-8 


-0.01433745 


1.167029 


X 


10" 


■10 


160 


1.009508 


-0.006329388 


-4.27 


X 


10" 


-8 


-0.01253953 


5.998161 


X 


10" 


■ 11 


180 


1.008438 


-0.005618167 


-2.66 


X 


10" 


-8 


-0.01114230 


3.334080 


X 


10" 


■ 11 


200 


1.007585 


-0.005050642 


-1.74 


X 


10" 


-8 


-0.01002523 


1.971428 


X 


10" 


11 


250 


1.006054 


-0.004032327 


-7.13 


X 


10" 


-9 


-0.008016120 


6.476438 


X 


10" 


■12 


300 


1.005037 


-0.003355744 


-3.43 


X 


10" 


-9 


-0.006677846 


2.607338 


X 


10" 


■12 



TABLE I: Summary of the gravitational self-force effects lU*/^, it"u''/ij^j and dE/dt for a variety of radii, Rn- The label — 3PN 
denotes the remainder after the post-Newtonian expansion in Eq. (2) has been removed from im*. 
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orbit about a large black hole of mass m. We use a 
standard frequency domain formalism that has been in 
continuous use for nearly forty years. This formalism is 
based upon a decomposition of the components of in 
terms of tensor harmonics. And we rely heavily upon the 
results of Regge and Wheeler [19] and of ZeriUi [14, 15]. 
We use the Regge- Wheeler gauge for h'^^ , although all of 
our results are independent of the gauge choice. 



for the regular field. 

We show in Sec. Ill that our problem only requires 
the regularization of u'^vf'hab- Following the original pre- 
scription of Barack and Ori [7] and extending it as in 
Ref. [23], we first evaluate 



ret(i) 
b 



I 

E 



Tct{lr, 
b 



(44) 



-I 



B. Regularization 

Two different approaches to regularization have been 
developed for self-force problems involving point particles 
and gravitational fields and both lead to the same con- 
clusion. Mino, Sasaki and Tanaka [3] follow the DeWitt- 
Brehme [2] analysis for electromagnetic fields in curved 
spacetime and show that the gravitational self-force may 
be described completely as geodesic motion in the per- 
turbed geometry gab + where is the part of 
the retarded metric perturbation which comes from the 
support of the Green function within the past null cone 
of the particle. Quinn and Wald [4] invoke a physically 
appealing "comparison axiom" and arrive at essentially 
the same conclusion. 

A third approach [5] notes that the tail part of the 
metric perturbation is only a portion of a solution of the 
perturbed Einstein equations. However a Green func- 
tion, different from the retarded Green function, deter- 
mines the singular part of the metric perturbation /i^^ 
which exerts no force on the particle. The remainder of 
the metric perturbation h^^ = — /i^^ is regular at the 
particle and is a source- free solution of the perturbed Ein- 
stein equations. The self-force is subsequently described 
as geodesic motion in the combined metric gab + h^i,. At 
first perturbative order gab + ^ vacuum solution 

of the Einstein equations, which has the desirable impli- 
cation that a local observer would see the particle move 
along a geodesic of a vacuum solution of the Einstein 
equations and, in fact, would only observe the combined 
field gab + h^b have no local method for distinguishing 
hfb from gab- 

A practical method for applying the regularization for- 
malism was described first by Barack and Ori [7-10] , later 
by Mino, Nakano and Sasaki [11, 12] and subsequently 
extended by others [20-25]. In this procedure the mul- 
tipole moments of /i^j and its derivatives are calculated 
analytically and referred to as regularization parameters. 
The sum of these moments diverges when evaluated at 
the particle, but each individual moment is finite. By 
construction /i^j, precisely matches the singularity struc- 
ture of the retarded field at the particle. Thus the differ- 
ence of the retarded moments and the singular moments 
gives a multipole decomposition of /i^^ which converges 
when summed over / and m. Schematically, this proce- 
dure gives 



ah 



R{lm) 



E 



ret(/m) 
b 



S{lm) 
^ab 



(43) 



Then we use the ansatz that 

C 



D 



~ ^ 1 + 1/2 [2l-l)[2l + i) 



El 



(2Z-3)(2Z-l)(2Z + 3)(2Z + 5) 
0{l-') (45) 



where B, C, D, Ei, . . . are the regularization parameters. 
The particular I dependence of the D, and Ei, . . . terms 
is related to the expansion of (1 — cos6')"+^/^ in terms of 
Legendre polynomials Pi{cos6); details are derived and 
described in Appendix D of Ref. [23] . 

For our problem it is known analytically that C — 
and 



y r'^yr — 2m) 



(46) 



where 2F1 is a hypcrgcomctric function, and r is the 
Schwarzschild radial coordinate of the circular orbit. 
Work in progress will describe the derivation of a va- 
riety of regularization parameters, including these. This 
knowledge of B and C, but not D, implies that 



ah 



rot(/) 
ab 



U U 



buS(l). 
"■ab ) 



(47) 



converges as 1/L To increase the rate of convergence, 
we augment our knowledge of B and C by numerically 
determining further regularization parameters [23]: We 



use the fact that the behavior of u°-u h 



b, ret(0 



evaluated at 



/i, must match as given in Eq. (45) for large 

This allows us to fit the numerical data to determine the 
additional regularization parameters D, and Ei to £'3. 
This results in the sum converging as ~ We sum up 
to ^ = 40, providing sufficient accuracy for the results as 
presented. 

A second gauge independent quantity is dE/dt as given 
in Eq. (19). The regularization parameters for this quan- 
tity are all zero, no regularization is required, and the 
sum over I converges faster than any power of I. 

Table I has a complete set of the interesting, 
gauge-independent data regarding a circular orbit of a 
Schwarzschild black hole. Any other gauge independent 
quantity, known by us, may be derived from these data. 
The results for iw* with the post-Newtonian expansion 
are displayed in Figs. 1 and 2 of Sec. I. 
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Barack and Sago [26] have recently used a time- 
domain formalism to calculate the actual gravitational 
self-force in the Lorenz gauge. Their Fr corresponds to 
our ^^dr{u°-u''h^i^y However, we worked in the Regge- 
Wheeler gauge and are not able to compare results. 



C. Gauge difficulties witli regularization 
parameters 

If the singular and retarded fields in Eq. (43) are in 
different gauges then the schematic description of regu- 
larization fails. 

The singular field is commonly described in the Lorenz 
gauge iz^ah' while the retarded field is easily found in 
the Regge- Wheeler [19] gauge rw^ab* ' ^^^'^ context of 
the Schwarzschild geometry. By gauge choice rwht^, = 
for even-parity perturbations, while iz/it^ is generally not 



zero and is singular at /i. In this case 



H41 



necessarily singular, and the "regularization" procedure 
fails. 

We now show that this failure is circumvented when a 
gauge independent quantity is calculated. 

Assume that the gauge vector relating the two gauge 
choices is known for the retarded field. 



(48) 



and is used to change the retarded field into the Lorenz 
gauge mode by mode, 



S{lm) 



iVGt{lm) iS{lm) 
Izh-ab - I'^'^ab ■ 



(49) 



Each component of izh^b now regular and the sum over 
modes converges. 

In this paper, we are eager to calculate the gauge in- 
dependent combination u'^u^h^jj for quasi-circular orbits 
which involves the sum over multipole terms such as 



R{lm) 



-.ar-.b ,rct(iji 
I U. ivjl^ab 



-a-b uS(lr> 



(50) 



The demonstration that u°'u^hab is gauge independent 
for quasi-circular orbits does not rely on hab being the 
retarded, the singular or the regular field, and it clearly 
shows that iZ"'«^V(a^^''"^ = at the particle, for any 

gauge vector^. Thus the Va^b term in Eq. (50) is identi- 
cally zero, and 



u^u^h^b 



r ret(^m) 



S{lm) 
b 



(51) 



^ The helical symmetry need not bo respected by , but the entire 
discussion of the gauge independence is then substantially more 
complicated. 



converges to a value that is independent of any gauge 
choice. 

This technique can be used to regularize any gauge- 
independent linear combination of components of hab and 
its derivatives. In fact, Moncrief [27] noted that the non- 
zero components of rw^ah may be described in terms of 
gauge independent, linear combinations of components 
of hab and its derivatives in a generic gauge. Thus, if 
the regularization parameters for rw^ah known and 
if /ijj^* is known in an arbitrary gauge, then what might 
be termed the regularized gravitational self-force in the 
Regge- Wheeler gauge could be determined in just this 
manner. And the result would, at least, be mathemati- 
cally well defined. 



V. CONCLUSION 

In Ref. [28] Barack and Ori state "The meaningful de- 
scription of the gravitational self-force must include both 
-^scif ^^'^ the metric perturbation hafs" We agree whole- 
heartedly with this statement. But we go a step further 
and believe that the value in calculating the gravitational 
self-force, in any particular gauge, is to apply it to a ques- 
tion whose answer is related to some physical observable. 
And a physical observable ought to be independent of the 
gauge choice. In Appendix C we give an example of how 
easily gauge confusion appears even in Newtonian physics 
where the direction of the gravitational self-force contains 
no useful physical information without additional knowl- 
edge about the coordinates and how the self-force is being 
evaluated. 

Similarly, we prefer to describe the effects of the self- 
force as geodesic motion in the perturbed and regularized 
metric gab + hab- The alternative description "acceler- 
ation agp in gab' depends upon the gauge choice and 
bears no relationship to any actual acceleration which an 
observer moving with the particle would measure with a 
collection of small masses and springs. At this level of ap- 
proximation fi is in free-fall in the actual, physical space- 
time metric gab + hab- To give Ogp in some gauge seems 
to imply that Ogp contains a true physical consequence — 
if that is the case we would prefer a description of that 
consequence. 

The circular orbits of the Schwarzschild metric perhaps 
provide the simplest framework for any self-force prob- 
lem. The angular decomposition and elementary discrete 
frequency spectrum imply that only ordinary differen- 
tial equations need to be solved in order to determine 
the metric perturbations. Apparently, few problems can 
be formulated in a gauge independent way within this 
simple framework. We would be surprised if another 
first-order self-force effect in this particular context were 
found which did not have a solution in terms of the data 
available in Table I. 

A general orbit has none of the natural symmetry of 
the quasi-circular orbit. Then none of our fi, u*, u"^, i? or 
J can be described in terms of gauge independent quanti- 
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ties. In this situation, where understanding the meaning 
of a question is likely to be as difficult as determining its 
answer, perhaps the only significant questions concern 
the waveform. 

However, our preference for Gauge invariance is a mat- 
ter of taste. So it is important to note that the value 



of 



'SF 



in the Lorenz gauge for circular orbits is pre- 



cisely defined and can be determined in a precise way, 
as ably demonstrated by Barack and Sago [26]. Their 
time-domain implementation is sharply focused on Ogp 
and takes a significant step in the direction of calculat- 
ing the elusive self-force effect upon waveforms. 

There is reason for optimism regarding waveforms. 
The difficulty of two disparate length scales for the 
extreme-mass-ratio problem can be avoided by using 
the analytically known singular field h^i^ to construct a 
smoothed out source [26, 29-31] for the wave equation of 
an approximation to the regular field The smoothed 
out source shows no structure with a length scale of fi. 
The numerically determined /i^j, can be guaranteed to be 
at least at /x [13, 31]. In this case h^^ directly gives 
the self-force effects on the motion and simultaneously 
provides the waveform in the wave zone. 

To take full advantage of the self-force formalism will 
require that /ij^^ be evaluated at second-order in ii. But, 
this impediment does not appear to be fundamental. It is 
likely that techniques based upon the current successes of 
numerical relativity will be able finally to reveal a better 
description of the dissipative effects of the self-force on 
gravitational waveforms. 



APPENDIX A: CONVENIENT IDENTITIES FOR 
GEODESICS OF THE PERTURBED 
SCHWARZSCHILD GEOMETRY 

The following identities arc useful for simplifying an 
equation by transforming terms involving E and J back 
and forth into terms only involving u°'. These equations 
are all elementary consequences of Eq. (9): 



u^-u^hab 



rE J 

-v^htb + -^u''h4.b, (Al) 



r — 2m 



du" 2mE 2J_, , , 

-TT—uhab = -7 7: — zT^u htb :ru h^b, (A2) 

or (r — 2m.Y r-^ 



and 



du°- 2(r-3m)^_^ 
or 



(r — 2m)^ 



Eu^htb, (A3) 



2m 



u''u''hab + 



r(r — 2m) 



dr 



2(r-3m) , 

-JuhA 



r^[r — 2m) 



(A4) 



APPENDIX B: GAUGE TRANSFORMATIONS 



The change in the metric perturbation under a gauge 
transformation .tJJ^,^ = x^j^j -f where ^" is considered 
to be infinitesimal (also known as an infinitesimal coor- 
dinate transformation) is 
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A/lafc = -C^gab = -S.^'dcgab - gacdbC " gcbdaS." 

= "^a^b-VbCa- (Bl) 

To preserve reflection symmetry we assume that and 
its derivatives are all zero on the equatorial plane. Then, 
for the Schwarzschild geometry we obtain 



2m 



(r- 277i)2^ (l-2TO/r) 



Ahtr = -drS,t 



2m 



re 



r{r ~ 2m) l~2m/r 



drC, (B2) 



dtC, (B3) 



2?77 



Aht4> = -dt^4, - 84,^1, 



and 



A/i, 



r 1 — 2m/r 



(B4) 



(B5) 



(B6) 



(B7) 
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For quasi-circular orbits of perturbed Schwarzschild we 
consider htt + ^h^t and ht^ + flh^^] these are propor- 
tional to u'^hat and W^hatj, respectively. Under a gauge 
transformation 

2m 

Ahtt + nAh^t = ~2ds -m^^-nd^(t, (B8) 

Combining these first provides 



2{dt + nd^)Ct - 2n{dt + nd^)C^, (bio) 



and subsequently 



U^U^Ahab = 



(1 - 2m/r)^ 



A{hu + 2nht^ + n'^h^^) 



(1 - 2m/r)2 



.(Bll) 



With the assumption that respects the helical sym- 
metry £/c^" = (dt + ild^)^^ = 0, the above effect of a 
general gauge transformation simplifies to 



U^U^Ahn 



E^ 



(1 - 2m/r)2 
0, 



2r 



(B12) 



when evaluated at the particle where 17^ ~ ra/r^ + 0{h). 
Equation (B12) and Am" = 0{h) then imply that 
A[u''vl'hab) = 0, at Oili). The fact that 



A {u'^u^'dthab) = 0, 



(B13) 



at the particle, follows simply from Eq. (Bll). A similar 
argument reveals the gauge independence of u'^u^d^hab 
when evaluated at the particle. 

It also follows from Eq. (Bll) that 



U^'vl' — Ahab = 

or 



E' 



(1 - 2m/ry 



—6m 



r (B14) 



evaluated at the particle. For a quasi-circular orbit 
£'2/(l-2m/r)2 = l/(l-3m/r) + 0(^), so that a general 
gauge transformation induces a change at the particle 



A{u''u''drhab) = - 



6m 



r'^{r — 3m) 



(B15) 



This result shows that the expression for fi^ in Eq. (28) 
is invariant under a general gauge transformation, while 
the radius of the orbit changes with the radial coordinate, 
''new = '"old + S,^ , and is not gauge independent. 



APPENDIX C: GAUGE DEPENDENCE OF THE 
NEWTONIAN GRAVITATIONAL SELF-FORCE 

In Newtonian gravity, when one mass mi of a circular 
binary is infinitesimal, then the orbital frequency is Vl^ — 
Gm2/r^. If mi is small but finite, then Q, changes by 
0(mi/m2). This change in 17 is properly described as a 
consequence of a Newtonian gravitational self- force [32] . 

The orbital acceleration of mi in a binary is given by 



Gmo 



in + r2)2 



(CI) 



where ri and r2 are the radii of the orbits and arc also the 
distances from the masses to the center of mass. Thus 
TOiTi = m2r2. 
Note that 



lim ri + r2 = ri . 

mi— *0 



And in this same limit, 



lim rifl'^ — Gm2/rl. 



(C2) 



(C3) 



But how does the acceleration of mi change as the limit 
is taken? The answer would determine the sign of the 
self-force. For an extreme mass ratio, the magnitude of 
the acceleration of mi is precisely 



on the one hand, but also 



Gm2 



(ri+r2)2 



Gm2 



Gm2 



{ri+r2y 



2m 1 
m2 



(C4) 



(C5) 



on the other. 

The acceleration and the self-force depend upon the 
detail of how ri + r2 approaches ri . If the distance ri + 
r2 between the masses is held fixed during the limit (ri 
grows slightly while r2 shrinks) then Eq. (C4) implies 
that the acceleration of mi is constant and that there is 
no self-force. However, if mi's orbital radius ri is held 
fixed during the limit, then Eq. (C5) implies that the 
limit of the acceleration is approached from below, and 
the self-force points outward. 

Even in Newtonian physics the gravitational self-force 
is ambiguously defined and not particularly relevant to 
understanding Eq. (CI). 

This Newtonian ambiguity is precisely equivalent to 
the gauge ambiguity of the gravitational self-force in gen- 
eral relativity. For an extremc-mass-ratio binary, if the 
origin of the Schwarzschild coordinates is at the center 
of the black hole then the coordinate value of r at the 
orbit would represent the distance between the particle 
center and the black hole center, i.e. ri + r2, and Eq. 
(C4) would imply that there is no self-force. However, a 
dipole gauge transformation can move the origin of the 
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coordinates to the center of mass of the system and then 
the coordinate value of r at the orbit would represent the 
distance between the particle and the center of mass, i.e. 
ri alone, and Eq. (C5) would imply that the self-force 
points outward. Consequently, a simple dipole gauge 
transformation substantially changes the appearance of 
the gravitational self-force. 

The resolution of this confusion in Newtonian physics 
is to refrain from introducing the concept of the self-force. 
It seems reasonable that general relativity would follow 
the Newtonian lead. 



APPENDIX D: USING THE ENERGY FLUX FOR 
RADIATION REACTION 



For the special case of quasi-circular inspiral, it has 
long been common knowledge that the dissipative effects 
of radiation reaction could be calculated by determining 
(i) the total energy of particle fiE, (ii) the orbital fre- 
quency n, and (in) the rate of energy loss dEtot/dt via 
gravitational waves. Each of these is to be calculated as 
a function of the radius r of the orbit of ^. As energy is 
lost fj, slowly spirals inward, and the rate of change of f2 
is then expected to be 



dn dEtot dQ /dr 
dt dt udE/dr' 



(Dl) 



In this Appendix we show that dfl /dt in this equation is 
equivalent to dil/dt as determined in Eq. (1), which is 
based upon the self-force formalism. 

Equation (Dl) and Eq. (1) appear similar but have an 
important difference. In Eq. (Dl), dEtot/dt is the rate 
that energy is radiated out at a large distance and down 
into the black hole, while in Eq. (1) jj. dE/dt is determined 
locally from hat via Eq. (19) and is the rate at which the 
the self-force removes energy from the particle. These 
two equations give the same dVl/dt only if the rate that 
energy is lost through gravitational waves is equal to the 
rate that the self- force removes energy from the orbit. 
We now show that this is the case. 

We assume that gab is the Schwarzschild metric, and 
only in this appendix an overdot represents a derivative 
with respect to the Schwarzschild time coordinate t. 

The perturbed Einstein equations, with a perturbing 
stress-energy tensor Tab = 0(/i) being the source, may 
be written as (see Eq. (35.58) of Ref. [33]) 



l&TlTab = -Eabih) 



(D2) 



where 



Eabih) = V'^hab + VaVbh-2Vi^a^''Ky 

+2Ra'b''Kd + gabiy'V'Kd - V2/i),(D3) 

with h = habg"^''- 

For arbitrary symmetric tensors fc°^ and hab, the op- 
erator Eab{h) satisfies the identity [13] 

k'^'^Eabih) = VcF^k, h) - hab) (D4) 



with 

F''{k,K) = k''''\/''hab ~ ^kVh - 2k''''V''hab (D5) 
where hab = hab ~ ■^Qabh'^cj and similarly for kab- Also 

- 2^ ak'^'^'^hbc - 2k''''Ra\''Kd (D6) 

is symmetric under interchange of hab and kab- 
Let kab = hab- Then a consequence of Eq. (D2) is 

W^T'^^hab = -VcF^h, h) + hab) (D7) 

We now evaluate an integral of both sides of this equation 
over a constant t surface E. 

The stress-energy tensor for a point mass moving along 
a path through space X'^t) is 



■.5'^{x' - X\t)). 



'9 

So, an integral of the left hand side of Eq. (D7) is 



(D8) 



16n / T'^'habT^ sine drd0d(f) 
Is 



where the right hand side is evaluated at the particle. 

Assume that S is bounded by one two-sphere SEqo at 
a large radius and by a second two-sphere 9E2m close to 
the event horizon. And integrate Eq. (D7) over S. Then 
substitute Eq. (D9) into the left hand side and the result 
is 

u'^u^h J f 
l6TTfi / y ^. F^ih.hy sin edrdOdcj) 



+ {h''\haby sine drdOi 



(DIO) 

where the time derivative was moved outside the second 
integral by the virtues of the time independence of the 
metric and of the symmetry of ( , ) . That same term is 
then zero because ^ is moving along a circular geodesic 
and hab respects the helical symmetry Ckhab = 0. Thus, 
the integral is 



167r/i 



u°-u''hab 

d 

~ ~di 



dc[r^ sin OF'' {h,h)] dr dO d<j) 
F\h, h) sin edrded(j) 
J -^[r'^ sin eF\h,h)] dr dO dcj) 
r^ sin eF^'ih, h)] de dcj) 



(Dll) 
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where the first term on the right hand side of the second 
equation is zero because of the hehcal symmetry, and the 
third equaUty follows from Gauss' law. 

The definition of F^ik, h) in Eq. (D5) gives 



h) 



h'^'^Vhab - - W/i ~ 2h''''W''hab. (D12) 



Note that when the integrals in the last line of Eq. (Dll) 
are evaluated close to the event horizon and far into the 
wave zone F^{h,h) is then equal to 327r7^'', a compo- 
nent of the effective stress-energy tensor of a gravitational 
wave in a background geometry as given in Eq. (35.70) 
of Ref. [33]. 

Thus, the right hand side of Eq. (Dll) is 327r times 
dEtot/dt, the rate that the gravitational waves carry en- 
ergy out through the boundaries of S. Our conclusion 
is 



dEt. 



dt 



U^U^dthab 



dE 



(D13) 



so that Eq. (Dl) and Eq. (1) have the same implications. 

In a numerical study of quasi-circular orbits, a com- 
parison of dEtot/dt with ^dE/dt is a useful test of the 
numerical implementation and accuracy, but is not a test 
of the self-force formalism. 

For completeness, we note that the actual, slow inspi- 
ral of a quasi-circular orbit is solely a result of the ra- 
diative energy and angular momentum losses which lead 
to n slowly falling inward to a tighter orbit. The radial 
component of the self-force does no work on fi and is not 
responsible for the inspiral. 



APPENDIX E: THE POST-NEWTONIAN 
EXPANSIONS OF u* AND Q 

Recently we have used the analysis by Blanchet, 
Faye and Ponsot (BFP) [17, 18] to determine the post- 
Newtonian relationship between the two gauge indepen- 
dent quantities Q and u* in the extreme mass ratio limit. 

We use many of the results and much of the notation 
of BFP without extensive clarification, and we limit our 
interest to orbits which are circular up to the effects of 
radiation reaction. In this appendix only, c is not set 
to unity, and an expression of the form 0{n) refers to 
general terms of order 1/c" in the limit c ^ oo. 

The two components of the binary system are mi and 
7712, and 



M = mi + TO2 



(El) 



is the total mass. The components are located at yi 
and 2/2 with velocities vi = dyi/dt and V2 = dy2/dt, 
whose magnitudes vi and V2 are given in Eqs. (E15) and 
(E16) below. For convenience the Cartesian dot product 
between two vectors is denoted 



{V1V2) = Vi-V2= -Vl V2 



(E2) 



for the circular orbits of interest. 

Also the relative position of the masses is ri2 = j7i — 2/2, 
and the relative velocity is V12 = vi—V2- The magnitude 
of the relative velocity is 



V12 = + 0(6) 



(E3) 



where f2 is the angular velocity. The 0(6) term is a 
consequence of radiation reaction at 2.5PN. 

We summarize some results from BFP and Blanchet 's 
Living Review [18]. One dimensionless measure of the 
strength of the gravitational field is 



GM 
ri2c2 



and a second is 



^12 



-0(8). 



(E4) 



(E5) 



These two parameters are precisely equal only in the 
Newtonian limit. 



7=^+0(4). 



(E6) 



It is natural, then, to introduce a third dimensionless 
parameter 



„3 _ „,2 "12 



a; = 7 — = 7 2 
c c 



(E7) 

which has the useful features that x = i'i2/c^ + 0(4) and 
that it depends only upon the quantities M and 51 which 
arc independent of any choice of coordinates. 

The second dimensionless parameter Wi2/c^ is known 
in terms of 7 at the 2PN order to be 

^ = 7 - (3 - y)!"^ + + + 1/^)7^ + 0(8) (E8) 

as given in BFP Eq. (8.6) or Eq. (190) in Ref. [18], with 
V = mimz/M^. From Eqs. (E7) and (E8) 

= 73 _ (3 „ j,)^4 _^ ^ ^ j,2^^5 ^ (^(^2). (E9) 

Alternatively, 7 may be expanded in terms of x which 
results in 

7 = X + (1 - .2 + (1 - ^) .3 + 0(8) (ElO) 

as in Eq. (193) of ref. [18]. 

Now we focus on m* and its post-Newtonian expansion 
via the regularized metric at mi given by BFP in their 
Eq. (7.6). 

Let V"' = (c, iTi), and let be parallel to but 
normalized so that u°-Ua = —1. Thus, = u^V"" is the 
appropriately parameterized four velocity of mi. In this 
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Appendix our major task is to find the post-Newtonian 
expression for 

u* = l/V-V'^ybgab = l/^J\gtt + 2v\g^t + vlv{g,j\. 

Our restriction to circular orbits allows for simplifi- 
cations of the more general regularized 2PN metric at 



the location of mi given in BFP Eq. (7.6). For exam- 
ple, (ni2i;i2) = ri2 ■ Wi2/fi2 = 0{5), and such terms in 
Eq. (7.6) do not lead to a contribution to u* at 0(6). 
Starting with BFP Eq. (7.6), specializing to the circular 
orbit case, and dropping terms whose contributions to u* 
are smaller than 0(6), we obtain the metric components 



(5oo)i = -1 + 2 



Gm2 Gm2 



Avi 



ri2 



- 2- 



ri2 



Gm2 



G^mims /23 2 39, 

^-i-y(^'i«2) 



12 



c^ri2 



ri2 



' 12 



' 12 



i6^2 



12 



47G^mim2 , Gm2 f Gm2 r on G^toi^ 17G^TOim2 G^TO2^ 

-Wo + — — < [-w^J ^ ^- ^ ^ 2 — ^ — 



' 12 



0(8), 



(E12) 



(5oi)i 



A ^7712 , 

-4^ V2 

c^ri2 



Gm2 
c^ri2 



Gtoi j 

4 vl 

ri2 



Gnii Gm2 



ri2 



ri2 



0{7) 



(E13) 



(5y)i 



&rx2 



Gm2 

^^2 



Gnii Gm2 



ri2 



ri2 



Gm2 



[AvW2}+0{6) . 



(E14) 



The order of terms, bracketing and other details in these 
equations are as close as possible to the original form 
given in BFP to facilitate a comparison of this form to 
the original, which an enthusiastic reader might attempt. 
In the center of mass frame of reference 

Vl =\vi\ = [7712 + 3-f^iy{mi - 7122)] V12/M + 0(6), (E15) 



and 

V2 = \v2\ — [mi - 37^1^(7711 — 7772)] V12/M + 0(6), (E16) 

from Eq. (187) of Ref. [18]. 

With these substitutions, the use of Eqs. (E12)-(E14) 
and G being removed in favor of 7 via Eq. (E4), we obtain 



i/u --V V gab - I- -TPnr - 2^rr - ^ ..^o \} + ^/^i + ^Mi) + ,.0 + ■5/^1) - 4- 



4 



(4 + 27^.1 + 114^2 ^ 47^^3) + !|V _ ^^^^ ^ 3^2) ^ ^(3)^ (g^^) 



with 111 z= 7771 /M. 

Next 7'i2/c2 is removed in favor of 7 using Eq. (E8), 



1/u*^ = 1 + (-3 + 4/ii - /i^) 7 + (3 - 10/ii + - 5/i? + /4) 7^ 



2 

And 7 is removed in favor of x using Eq. (ElO) which yields 



1^1 + - |m? - f + 6m? - M?) 7^ + 0(8). (E18) 



1 , 3 „ 1 2 \ / 27 13 53 2 4 . 1 4 , 



/135 175 409 97 69 1 1,6^3,0^8^ 



as the 0(6) expression for u* of r77i in a circular binary with 7772. 
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To change from this expansion, appropriate for a com- 
parable mass binary, to a similar one with an extreme 
mass ratio it is necessary to replace the post-Newtonian 
dimensionless parameter 



Q.Gm2 I ^ mi 

1712 



with a dimensionless parameter 



2/3 



(E20) 



(E21) 



more convenient for simultaneous, independent expan- 
sions in v/c and in e = TOi/to2- [Note that = 
Gm2/Rac^. from Eq. (36).] It follows that 



a; = Xe(l-He)2/3, 
and expanding yields 



2 1 , 4 , 

+81^ 



(E22) 



(E23) 



The final steps are the substitutions of /ii = e/(l + e) and 
of Eq. (E23) into Eq. (E19) and the subsequent collec- 
tion of powers of the small quantities e and of . Thus, 
for an extreme mass ratio with a particle of mass mi 
in a circular orbit with the orbital frequency f2 about a 



Schwarzschild black hole of mass m2 



28 
27^ 



135 



97 , 97 



5e + — - —e-" x\ 



16 8 6 

+0(8) + 0(2) X 0{e^). 



3 \ „3 



(E24) 



This result is consistent with the post-Newtonian expan- 
sion of the geodesic value of u* 

ou* = (1 - 3m/r)-i/2 = (1 - ix,)-^''' 

= l + ^., + |x2^^^3_^0(.,)4. (E25) 

And the first order self-force effect on u* is 

lu* = ~ex, - 2exl - hex^ + O(ex^), (E26) 

which is plotted in Figs. (1) and (2). In the future, a 
second-order perturbation analysis might be compared 
with 

2u' = ^x, + ie^xl+^-l^xl + 0{e^xt). (E27) 
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